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Abstract 

The normal muon decay fi + — > e + v e T7^ is studied as a tool to discriminate between the Dirac 
and Majorana type neutrinos. For the purpose to do this, we propose a new parameterization in 
place of the Michel parameters for the energy spectrum of e + . The x 2_ fitting is used by noticing 
different energy spectra between the Dirac and Majorana neutrino cases. We assume the interaction 
Hamiltonian which consists of the V — A and V + A charged currents. 
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I. INTRODUCTION 

Neutrinoless double beta decay [l| is a best subject for investigating whether the neutrino 
is of the Dirac or Majorana type. However, it needs to improve the present sensitivity of the 
detector for observing this decay mode, because of the tiny neutrino masses and/or the small 
contributions from the right-handed weak current. Under this circumstances, it seems to be 
meaningful as a complementary study to survey whether muon decay can be used as a tool 
to determine the type of the neutrino. Recently, we have proposed a new parameterization 
of the e ± energy spectrum of the muon decay that is suitable for discriminating between the 
types of neutrino 2]. It is shown that there is a possibility to take advantage of the different 
energy dependences for their energy spectra. We propose a method in which the \ 2 value 
determined experimentally by assuming the Dirac-type neutrino is compared with the one 
determined for the Majorana- type neutrino. This may provide a test to determine the type 
of neutrino, although it is indirect. 

II. GENERAL FRAMEWORK 

We assume the following effective weak interaction Hamiltonian for the ^ decay P 

H w (x) = ^{jl LaJ « L + XjI r Xr + vArJZl + kJIlJSr} + H - c - > (2-i) 

where Gf is the Fermi coupling constant. This weak interaction is naturally expected from 
the gauge models that contain both V — A and V + A currents with the left- and right- 
handed weak gauge bosons. The appearance of the coupling constant A is mainly due to 
the right-handed weak gauge bosons Wr, while terms with rj and k come from the possible 
mixing between the left-handed and right-handed weak gauge bosons, Wl and Wr. In the 
SU{2) L x SU{2) R x U{1) gauge model, the coupling constants n and rj become identical. 
However, they are treated as independent constants in this note in order to allow comparison 
with the more general case without a restriction from the gauge theory (see, e.g., |4(). 

The left-handed and right-handed charged weak leptonic currents, jn and jiR with flavor 
i = e and /i, are expressed in terms of the mass eigenstates of charged leptons Ee and 
neutrinos Nj with mass rrij as 

2n 2n 
JiLa = J^^lail ~ -YsWtjNj, jiRa = J^7a(l + l^VtjNj, (2.2) 
3=1 J=l 



for the case of the n generations. Here Uej and stand for the left-handed and right-handed 
lepton mixing matrices, respectively. 



III. DIFFERENTIAL DECAY RATE FOR NORMAL MUON DECAY 

The fj^ decay takes place as 

fj± _> e ± + Nj + W k , (3.1) 

where N k represents an antineutrino for the Dirac neutrino case, but it should be understood 
as N k for the Majorana neutrino case. 

If the radiative corrections are not included, the differential decay rate for emitted 
positron in the rest frame of polarized fi + is expressed as 

d 2 r(/i+ -> e+uV) f mfl G 2 F W 4 ' 



dxdcosO V 6-4 (V) 3 
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Jx 2 -xlD(x, 6), (3.2) 



where 



E m* a rn?, + writ 

x = Xo = -JL = 9.65 x lo- 3 , W = — ^ = 52.8 MeV. (3.3) 

W W 2m M 

Here m M and m e are the muon and electron masses, respectively, and E is the energy of e + . 
The angle 9 represents the direction of the emitted e + with respect to the muon polarization 
vector P M at the instant of the /z + decay. The allowed range of x is limited kinematically as 

x ^x^x max = (l-r 2 jk )~l with r%= (3-4) 

Here and m k are masses of neutrinos emitted in the /i + decay. 

The constant A in Eq. (|3.2j) is introduced to simplify the expression for the energy spec- 
trum by taking the arbitrariness of its normalization. It is referred to as a normalization 
factor. There are various possibilities for the choice of A, when experimental data are ana- 
lyzed, although these choices differ only by rearrangements of the terms in the theoretical 
expression. See Appendix A. 

IV. ENERGY SPECTRUM OF POSITRON 

We ignore some small terms proportional to both m e and neutrino masses (rrij and m&) 
in this note, for simplicity. Then, the e + energy spectrum part, xD(x, 6), is expressed as 



xD(x, 9) = xlN^ + P^ cos9P(x)}, (4.1) 

where P^ = \P^\ is the rate of muon polarization, and the isotropic part N(x) and anisotropic 
part P(x) are 

N{x) = - [a + {3x - 2x 2 ) + 12( k +c + e m k +m ) x (1 - x)] , (4.2) 
P(x) = \ [a^(-x + 2 x 2 ) + 12 (k_ c + e m k_ m )x {I -x)]. (4.3) 

Here, the decay formulae for the Dirac and Majorana neutrinos are obtained by setting 
e m = and e m = 1, respectively. The first terms with setting A = a + = a_ = 1 in these 
N(x) and P(x) correspond to the theoretical predictions from the standard model. The 
well-known Michel parameterization ^| is obtained if we choose the normalization factor 
A = Aio = a + + 2 k +c . 



V. COEFFICIENTS 

In the Dirac neutrino case, the coefficients in Eqs. ()4.2j) and ()4.3j) are expressed as follows: 

a ± =(l±A 2 ) and k± c = (k 2 ± rj 2 ) /2. (5.1) 

Here it is assumed that all neutrinos can be emitted in the muon decay. Then, the unitarity 
relations for the lepton mixing matrices U and V, namely, Y,j\Uej\ 2 = Ej|y^| 2 = 1 have been 
used. 

On the other hand, in the Majorana neutrino case, the coefficients are: 
a±= [(l-TT e 2 ) [l-u; 2 )±\ 2 - e 2 % 2 ], 

k± c = [K 2 (l-TT e 2 )% 2 ±r l 2 - e 2 (l-T^ 2 )] /2, (5.2) 
k± m = [k 2 \w^\ 2 ±r] 2 \w^J l \ 2 ] /2, 

where W 2 , W 2 , and w^Th are small quantities. Note that k± c has the same order of 
magnitude as k± m , in contrast to the Dirac neutrino case. In the Majorana neutrino case, 
we assume the existence of heavy Majorana neutrinos, which are not emitted in the muon 
decay. Then, we have Ej|[/^| 2 = 1 — W 2 and Sj|Vfj| 2 = U} 2 where the primed sums are 
taken over only the light neutrinos. In addition, the following products of U and V appear: 
= Sj U ej Vpj and = Y,' k V ek U^ k . 
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VI. ENERGY SPECTRUM OF e+ IN THE ji + DECAY 



We propose a new parameterization that directly represents deviations from the standard 
model. Namely, if we assume the SU(2)l x SU (2)#x U(l) model for simplicity, our expression 
for the energy spectrum of e + becoms 

xD(x, 9) = x 2 [(3 - 2x) + 2p c (3 - Ax) + \2e m p m {\ - x)\ 

+ x 2 P^cos9(-l + 2x), (6.1) 

where the parameters p c , p m , and £ are, respectively, given by 

Pc = -r— > 0, p m = — — > 0, and £ = , (6.2) 

^10 A-iQ Aiq 

with the choice of normalization factor Aio = a + + 2 k +c . Note that we have 



1 - A 2 
1 + X 2 + 2 V 2 

£ ~ 1 for the Majorana neutrino case, (6.4) 



£ = — for the Dirac neutrino case, (6.3) 



£ = 1 for the standard model. (6.5) 



A. Dirac type neutrino case 

Since the well-known Michel parameter pu is related to our p c as 2p c = (l — |pa/), we 
have the following expression: 

xD(x, 9) = x 2 [(3 - 2x) + 2p c (3 - Ax) + P M £ cos 9 (-1 + 2x)\ (6.6) 

= 6x 2 [(l -x) + Ipm{Ax - 3)] + P M £ cos#a; 2 (-l + 2x). (6.7) 
9 

The TWIST group 6] reported a precise experimental result 

p M = 0.75080 ± 0.00032(stat) ± 0.00097(syst) ± 0.00023, (6.8) 

by using angle-integrated spectrum derived from Eq.(6.7). The mean value is pu = 0.75080, 
although pm < 0.75 is satisfied within experimental uncertainty. Note that theoretical 
prediction from our model is pu = 0.75(1 — 2p c ) < 0.75, because p c = rf /(l + A 2 + 2rf ) > 0. 
They determined this pu by minimizing the \ 2 value. 
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B. Majorana type neutrino case 



The following expression is derived for the Majorana type neutrino case: 

xD(x, 9) = x 2 [(3 - 2x) + 2p c (3 - Ax) + 12p m (l - x) + cos 6(-l + 2x)\ . (6.9) 

Theoretically we predict 

1 1 

Pc ^ g ( k2 ^ 2 + ^ 2 ^ 2 ) > °. Pm - g ( k2 |«vI 2 + ^7 2 l«^| 2 ) > 0. (6.10) 

Therefore, there is a possibility to take advantage of the different x dependences of the 
terms including the parameters p c and p m in the energy spectrum, when y 2 value for the 
Dirac-type neutrino case is compared with the one for the Majorana-type neutrino case. 
This may provide a test to determine the type of neutrino, although it is indirect. 

VII. SUMMARY: HOW TO DETERMINE THE TYPE OF NEUTRINO 

There is a method that might make it possible to distinguish between two neutrino types. 
This method makes use of the different x dependences of the terms including the parameters 
p c and p m in the energy spectrum by comparing the x 2 values for the Dirac-type neutrino 
with those for the Majorana-type neutrino. Here we use the data with 9 = 7r/2, so that it 
may need new data-taking. 

For example, suppose we analyze experimental data with 9 = n/2 by using 

xD(x, 9 = tt/2) = x 2 [(3 - 2x) + 2p c (3 - Ax) + 12p m (l - x)} , (7.1) 

and obtain some x 2 value, say, Xm f° r the Majorana neutrino case. Then, suppose we repeat 
a similar analysis using 

xD(x, 9 = tt/2) = x 2 [(3 - 2x) + 2p c (3 - Ax)] , (7.2) 

and thereby determine Xd f° r the Dirac neutrino case. 

If Xm is much smaller than x\i we can conclude that there is a higher probability that 
neutrinos are of the Majorana type. 
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APPENDIX A: NORMALIZATION FACTOR A AND THE x 2 -FITTING 



We summarize the meaning of the freedom for the normalization factor A in the x 2 -fitting 
of the e + energy spectrum xD(x). 

First let us consider the case where xD{x) consists of the leading term f(x) and its small 
deviation term g(x) with weight w, namely, 

1 



xD(x) = — [f(x) + w g(x)\ 



(Al) 



Then the unknown constant w is determined experimentally by minimizing the x A -value 
that is, 

2 



a, 



Xa = s * 



d(x k ) - y A ,k 



o k 



(A2) 



where x k stands for the energy of the fc-th positron observed by experiment, and d(x k ) and a k 
mean, respectively, the corresponding experimental value of spectrum and its experimental 
error. The theoretical expression is y A>k = f{xk) + w g(xk) in this case. The summation is 
taken over all observed data. This case is referred to as analysis A, which corresponds to 
the choice of the normalization factor A = 1. 

Since we have \\ = a^w 2 — 2bpjv + ca in this analysis A, the minimal value of \\ anc ^ 
the corresponding value of w, respectively, are determined by 



Xa,, 



ca and wa = I — ) . 

a A J \a A/ 



(A3) 



Here a A , b A , and c A are defined by 

2 



ca = S fc 



Ok 

d( x k) - f(Xk) \ 

Ok 



b A = S fc 

2 



g(xk) \ ( d(x k ) - f(x k ) \ 



Ok 



o k 



(A4) 



Next, let us consider the case where terms in xD(x) are rearranged as follows, 



xD(x) = — [f(x) + wg(x)] 



1 + w 
A 



f(x) + 



w 



1 + w 



[g(x)-f(x)} 



(A5) 



Then we can determine ( j^) experimentally by minimizing the XB-value, 



2 v 



d{x k ) - y B , k 



Ok 



(A6) 
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where the theoretical expression is defined by 

w 

VB,k = f(x k ) + -^-^ \gM - f{x k )}. (A7) 

This case is referred to as analysis B, which corresponds to the choice of the normalization 
factor A = (1 + w). The minimal value of x| and the unknown constant w in this analysis 
B become as follows 

2 / „ °B 



XB, min =[cB-- , and ^=1—^). (A8) 



Here qb, 6b, and cb are defined by 



a B 



9{xk) ~ f{xk)\ 



J • 

In general, the value of u>b is different from wa- Therefore, the final value of w determined 
experimentally should be chosen as the wb, if Xa mm > Xb mini or v i ce versa. So, this means 
that there is freedom in the x 2- fitting such as the choice of the normalization factor A in 
the analysis of xD(x). 

It is worthwhile to note that the x dependence of [d(x) — f(x)] plays an important 
role in the above procedure. That is, in the muon decay, it means the x dependence of 
the difference between the experimental data and the prediction from the standard model. 
Let us consider the following imaginative example: If the x dependence of [d(x) — f(x)] 
happens to be proportional to that of g(x), that is, [d(x) — f(x)] = L^gix), then we have 
x\min = and correspondingly u>a = in analysis A. In other words, we have the relation, 
xXmin < xl,min- While, if [d(x) - f(x)} = L B [g(x) - f(x)}, then we have xl,mm = and 
correspondingly wb — L>b/(1 — Lb) in analysis B. Namely, we have the opposite conclusion, 

XA,min XB,min' 
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